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The proton polarizability correction to the Lamb shift of electronic and muonic 
hydrogen is calculated on the basis of isobar model and experimental data on the 
<D ' structure functions of deep inelastic lepton-nucleon scattering. The contributions 



of the Born terms, vector-meson exchanges and nucleon resonances are taken into 
account in the construction of the photoabsorption cross sections for transversely 
and longitudinally polarized virtual photons cry l- 
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5 ■ I- INTRODUCTION 

^P H| The precise investigation of the energy levels of hydrogenic atoms (muonium, positron- 

ium, hydrogen atom, deuterium, helium ions et al.) allows to obtain more exact values for 
many fundamental physical constants such as the lepton masses, the ratio of lepton and 
proton masses, the fine structure constant, the Rydberg constant which are used for cre- 
ating standards of units 0, 0- The insertion of new simple atomic systems in the range 
of experimental investigation can lead to significant progress in solving of these problems. 
Muonic hydrogen (/ip) is just one of a number of simple atoms which has attracted consid- 
erable attention in the last years. Since the muon is heavier than the electron by a factor of 
206, the energy structure of (/ip) is extremely sensitive to the effect of the electron vacuum 
polarization, recoil effect, the proton structure and polarizability corrections 0,0, E|- In 
the Lamb shift (2P — 2S) of (fip) the value of the proton structure correction of order (Za) A 
increases essentially as compared with electronic hydrogen. Considering that this contri- 
bution is determined also by the proton charge radius r p , the experimental investigation of 
the energy spectrum of muonic hydrogen can play the key role in a more precise study of 
the proton structure along with the experimental data on the electron-nucleon scattering 
@, 13 13- At present time the Lamb shift (2P — 25*) experiment in muonic hydrogen at PSI 
(Paul Sherrer Institute) with a precision of 30 ppm entered the closing stage (H, EJ. The 
suggested aim is to determine the root-mean-square (rms) charge radius of the proton to the 
accuracy 10~ 3 , about a factor 20 times better than presently known from electron scattering 
experiments. It opens the possibility to check bound state QED predictions toward a level 
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of 10 -7 precision. Because the Lamb shift experiment in (/u,p) measure the energy difference 
between the 2 3 S , 1 / 2 and 2 5 P^/^ atomic levels, it is important to know the hyperfme structure 
of 2S and 2P energy levels p,lU|- I n this connection it is important to consider again several 
theoretical contributions which have essential role in the calculation of the total theoretical 
value of the (2P — 2S) Lamb shift in the atom (u,p) with the neces sary accuracy. The proton 
polarizability contribution is among important effects 3, 111 3 U4lJJ,ll2l,llS|, which was 
calculated earlier in muonic hydrogen on the basis of experimental data on the structure 
functions of the lepton-nucleon scattering. Slightly different values for this contribution 
occur to the present 0, 17. lif. 

The study of electromagnetic excitations of baryonic resonances which is carried out at 
CEBAF (Continuous Electron Beam Accelerator Facility), entered currently a new phase of 
development because new data with unprecedented precision have become available. The 
aims of the CEBAF investigations are in the measurements of the nucleon transition form 
factors to nucleon resonances N* at different photon virtualities Q 2 in the resonance re- 
gion, in the study of the gluon content of baryonic resonances and the helicity amplitudes 
A1/2, A 3 /2, S1/2 for different states N* [lij . Recently, precise measurements for exclusive 
electroproduction of 7r°, 7r + , t] mesons on protons in the resonance region were obtained at 
Jefferson Lab using the CEBAF Large Acceptance Spectrometer (CLAS) The obtained 
experimental data include measurements of the cross sections, angular distributions for the 
7m + , pr) final states. These data allow also to investigate the contributions of resonances 
to 7r, r] electroproduction in detail. The performed experiments permitted already to refine 
the values of a number of theoretical parameters determining the production of the nucleon 
resonances N* in the reaction 7* + iV — > N*. The goal of the present investigation is to per- 
form new calculation of the proton polarizability contribution to the Lamb shift in hydrogen 
atom based on theoretical model (unitary isobar model), describing the processes of photo- 
and electroproduction of tt, rj mesons, nucleon resonances on the nucleon in the resonance 
region, and on evolution equations for the parton distributions in deep inelastic region. 



II. GENERAL FORMALISM 



The proton polarizability contribution to the Lamb shift of order (Za) 5 is determined 
by the amplitude of virtual Compton forward scattering 7* + p — > 7* + p, presented on the 
Feynman diagrams in Fig.l. It's parameterization has the following form |2lL \22\: 
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where k is the four-momentum of the virtual photon, v = k is the virtual photon energy, 
m 2 is the proton mass. Symmetrical part of the tensor (1) gives the contribution to the 
Lamb shift (the structure functions C\,i{y, k 2 )) and antisymmetric part contributes to the 
hyperfme structure (the structure functions Hi^iy, k 2 )). The lepton tensor is the following: 
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FIG. 1: Two-photon Feynman amplitudes determining the correction on the proton polarizability 
in the Lamb shift of hydrogen atom. 



where pi j2 , <?i,2 are four-momenta of the lepton and proton in the initial and final states, mi 
is the lepton mass. 

Taking the product of the amplitudes (1) and (2), we can extract the part of the inter- 
action operator which is not dependent on the particle spins and contributes to the Lamb 
shift (LS): 
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The structure functions Cj(/co, k ) satisfy a following dispersion relations [2 
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The threshold value of the photon energy u represents the minimal energy needed for the 
production of the 7r-meson in the reaction 7* + p — > tt° +p. Let us to point out that reliable 
data on the subtraction term in the first dispersion integral (4) are absent. But in the limit 
of small values of k 2 this term is connected with the proton magnetic polarizability: 
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where 0m — 1.9(0.5) x 10 4 fm 3 24]|. The dipole parameterization for /?m(^ 2 ) was suggested 
inRef.Q: 

A 8 

P M (k 2 ) = ^M (A2 + fc2)4 , (7) 

where A 2 = 0.71 GeV 2 as for the elastic nucleon form factor. Imaginary parts of the 
amplitudes Ci(ko,k 2 ) are expressed in terms of the structure functions Fi(x,Q 2 ) for deep 
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inelastic scattering as 

-ImCi(x,Q) = , -ImC 2 (x,Q) = , x = . (8) 

7T 777,2 V 2m 2 V 

So, to obtain numerical value of the proton polarizability correction we can use experimental 
data for the functions F^p iv, k 2 ) and different parameterizations for it prepared on the basis 
of these results IHSEj. Making use of relations (5)- (7) and transforming the integration 
in the loop amplitudes to four- dimensional Euclidean space with the aid of the formula 

d A k = Att / k 3 dk / sin 2 </> • d(p, k° = k cos 0, (9) 



we can perform the integration over the angle variable <fi and represent the proton polariz- 
ability contribution to the Lamb shift of hydrogen atom in the form: 

(fl+ w- t) {" 8V;(1+s)3/ ^ (2s+ii) " 

-Vi(t - As) [t + s(6 + 2R + As + t)} + VA + t[(t -2)t + s(8 + t 2 + 2R(t + 4))]|f 2 (t/, k 2 ), 

W = i + (i-|)(^-i), « = = m 

where R(y, k 2 ) = ol/^t is the ratio of the cross sections for the absorption of longitudinally 
and transversely polarized photons by hadrons. In such a way, the correction AE^ is 
expressed in terms of two structure functions F 2 (u, k 2 ) and -R(z/, k 2 ), describing unpolarized 
lepton-nucleon scattering. 



III. STRUCTURE FUNCTIONS OF UNPOLARIZED LEPTON-NUCLEON 

SCATTERING 



The greatest contribution to the integral (10) is given by the region of the variable k 2 : 
-h 1 GeV 2 and near the threshold values of a photon energy v. So, the exact construction 
of the structure functions F 2 , R in this region is very important to obtain reliable estimation 
of the proton polarizability effect. Deep inelastic lepton-nucleon scattering is described by 
the following reaction: 

l + N^l' + X, (13) 

where X denotes the sum over all final particles. This reaction represents the inclusive lepton 
(I) production with the measurement of their energy and scattering angle. It is assumed 
that this process occurs due to one-photon exchange. An important kinematical variable of 
the reaction (13) is the invariant mass of electroproduced hadronic system W: 



W 2 = m 2 2 -Q 2 + 2m 2 u, k 2 = -Q 2 . 



(14) 



FIG. 2: The proton polar izability correction in the resonance region. On the Feynman diagram 
solid, double solid, wave and dashed lines correspond to the nucleon, baryon resonance, photon 
and pion respectively. 



Using the variable (14) we can divide the total integration region in Eq.(lO) on the resonance 
region W < 2 GeV where the production of low-lying nucleon resonances occurs and deep 
inelastic region when W > 2 GeV. 

In the resonance region the proton polarizability contribution to the Lamb shift is 
determined by the processes of a 7T-, ?7-meson production on nucleons and the produc- 
tion of basic low-lying nucleon resonances. Several amplitudes of such reactions are pre- 
sented in Fig. 2. To calculate the contributions of separate resonances to the cross sec- 
tions orrin the isobar model we used the Breit-Wigner parameterization suggested in 
Refs. IMB ISiiSiiililiiHHB- In the considered region of the variables 
k 2 , W the most contribution is given by five resonances: P 33 (1232), S , n(1535), Pi 3 (1520), 
Pn(1440), Fi5(1680). Accounting the resonance decays to the Nn- and iV^-states we can 
express the absorption cross sections cri/ 2 and 03/2 as follows: 



^ 2 r 7 r R _ 
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where Ax/2,3/2 are transverse electromagnetic helicity amplitudes, 
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X = 0.3 GeV. 



The resonance parameters Tp, Mr, ji, J2, ^tot are taken from Refs.[2J, |3jj 
accordance with Refs. the parameterization of a one-pion decay width is 
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for the resonance P 33 (1232) and 



R 



r 1 (&±f 

1 \q 2 + 5 2 



(18) 



for P 13 (1520), Pn(1440), F 15 (1680). I is the pion angular momentum, 5 2 = (M R — m p — 
m 7r ) 2 +r|j/4. Here q (k) and q R (k R ) denote the c.m.s. pion (photon) momenta of resonances 
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with the mass M and Mn respectively. In the case of 5n(1535) we take into account ttN 
and r/N decay modes J32I Hsf : 



2 /"*2 

_ Q*,v u r <?7r >'? + ^, 

R-*Nn,N V ~ —°Tr,V L R o T r2 



w v:* (19) 



where is the tt, tj branching ratio. The cross section is determined by an expression 
similar to Eq.(15) where we must change Ax/2,3/2 on the longitudinal amplitude S\/ 2 - The 
calculation of helicity amplitudes Au^ A^/o , Sy 2 & s functions of Q 2 was done on the basis 
of the oscillator quark model in Refs. |37l. Eol 41, 0, H3, 3 • 



The two-pion decay modes of the higher nucleon resonances 5*11(1535), £>i3(1520), 
Pn(1440), Fi 5 (1680) are described phenomenologically using the two-step process as in 
Ref . |36j| . The high-lying nucleon resonance R can decay first into iV* (P33 or Pu) and a pion 
or into a nucleon and p, a meson. Then the new resonances decay into a nucleon and a pion 
or two pions: 

fiV* + 7r^iV + vr + 7r, 

\p(a) + N ^ N + n + n. v ' 

The total decay width of such processes can be presented as a phase-space weight integral 
over the mass distribution of the intermediate resonance r = N*, p, a (a = tt, N): 

r (W\-^L [ W ~ ma j 2 p 2 V r ,tot{p) (M R -m 2 -2m^) 2 + C 2 

L R ^ r+a [W)- w ^ a ^Pf 7r{fi 2_ m 2 ) 2 +lx 2 T 2 tM {W _ m2 _ 2m7r) 2 + C 2> W 

where C = 0.3 GeV. The factor P 27T must be taken from the constraint condition: 
r ' R-^r+a(Wjt) coincides with the experimental data in the resonance point, pf is the three- 
momentum of the resonance r in the rest frame of R. T ryto t is the total width of the resonance 
r. The decay width of the meson resonance in Eq.(21) is parameterized similarly to that of 
the P 33 (1232): 

(\ 2J r +l 2 1 r2 
i) ^ ' = < 22 > 

where m r and p are the mean mass and the actual mass of the meson resonance, q and q r 
are the pion three momenta in the rest frame of the resonance with masses p and m r . J r 
and r r are the spin and decay width of the resonance with the mass m r . 

Main nonresonant contribution to the cross sections <Jt,l i n the resonance region is de- 
termined by the Born terms constructed on the basis of Lagrangians of 'jNN, ^tttt , ttNN 
interactions. Another part of nonresonant background comprises the t- channel contribu- 
tions of p, to mesons obtained by means of effective Lagrangians / -yirV , VNN interactions 
(V = p,uj) (3^, Hij]. In the unitary isobar model accounting the Born terms, the vector 
meson, nucleon resonance contributions and the interference terms we calculated the cross 
sections <jt,l by means of numerical program MAID (http://www.kph-uni-mainz.de/MAID ) 
in the resonance region as the functions of two variables W and Q 2 . Then the structure 
function F 2 (W,Q 2 ) can be presented as follows: 

where K is the flux factor of virtual photons. There are two definitions for the quantity K: 
the Gilman definition Kg = VQ 2 + 1,1 1 an d the Hand definition Kh = ^(1 — Q 2 j1vn 2 v) = 
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FIG. 3: Total photoabsorption cross section crtot(W, Q 2 ) 
variables Q 2 (0 -r- 1) GeV 2 and W (1.1 -r- 1.8) GeF. 



{gt + gl) in as the function of 



(W 2 - m 2 )/2m 2 |4 
^(^,Q 2 ) = (a T +a L ) 



Our results are presented in Fig. 3 for the total cross section 
It contains three clear peaks corresponding to resonances P33(1232), 



.Di3(1520), Fi 5 (1680). The photoabsorption cross section at„t{W, Q = 0) differs from exper- 
imental data in the range 1.5 < W < 2 GeV obtained in Refs. |47L l48| . So, theoretical model 
must be improved by the account of two-pion resonance decays as described in Eq.(21). 

In the nonresonant region there exist several parameterizations for the function F 2 (Q 2 , W) 
(2o| obtained on the basis of experimental data on deep inelastic lepton-nucleon scattering. 
The structure function F 2 can be expressed in terms of parton distributions 



F 2 (x,Q 2 )=^e 2 xq t (x,Q 2 ). 



(24) 



So, to construct it in deep inelastic region we can use the Q 2 evolution equations for the 
quark and gluon distributions [49j | : 
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(25) 



(26) 



where the sum is considered over all quarks and antiquarks. P qq , P gq , P qg , P gg are the 
quark-gluon splitting functions Numerical solution of the integrodifferential evolution 
equations (25), (26) by means of the method suggested in Ref.|5l| allows to obtain the parton 
distributions and the structure function F 2 (x, Q 2 ) for different values of a photon momen- 
tum squared Q 2 . Corresponding numerical results are in good agreement with the world 
experimental data. On the other hand, the 23 parameter model for the structure function 
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F 2 (x,Q 2 ) based on experimental data was proposed in Ref. [2 71]. Here it was expressed as a 
sum of the Pomeron and the Reggeon F^- term contributions: 

Q 2 ) = 7^r^ Wi^ Q 2 ) + Q 2 )] , (27) 

t^/ ■+" TTIq 

F?(x,Q 2 ) = C n {t)x°£ {t \l - x) h ^\ F?(x,Q 2 ) = C P (t)x a / {t \l - x) b ^\ (28) 
where x is the Bjorken variable, 

1 W 2 -m 2 1 W 2 -m 2 

— = 1 + 7^— -r> — = i+ n2 ^ r » 29 



More accurate numerical values of the model parameters are presented in Ref.|27J (version 
2 (2004)). 

For the second structure function R(Q 2 , W) in nonresonant region we used also the pa- 
rameterization obtained with the aid of experimental data ([52]. In the resonance region 
there are no experimental data for the quantity R{Q 2 , W). In the most important part of 
the resonance region <jl <C ot- So, to perform the numerical calculation of the correction 
AE^ we supposed that R(Q 2 , W) « 0. 



IV. NUMERICAL RESULTS 



In this paper we calculate the proton polarizability correction to the Lamb shift of elec- 
tronic and muonic hydrogen on the basis of the isobar model describing the processes of 
low-energy scattering of virtual photons on nucleons and the evolution equations for the 
parton distributions. These two significant ingredients of the calculation allow to construct 
the absorption cross sections for transversely and longitudinally polarized photons by nu- 
cleons <tt,l and to express the structure function F 2 (Q 2 ,W) (23), (27), which determines 
required contribution (10). Numerical results are presented in Table 1. We investigated 
contributions to the correction (10) which have numerical value of order 1 fieV for the 15* 
state in muonic hydrogen and 1 Hz for the IS* state in electronic hydrogen. It is inferred 
from these results that the basic contribution to the polarizability effect is given by pro- 
cesses of the 7r-meson production on nucleons in the reaction 7* + N — > tt + iV including 
the resonant reactions (2 line in Table 1). In the isobar model we kept also processes of 
the ?7-meson production on nucleons [45| (3 line in Table 1) and the production processes 
of the .fT-mesons (4 line in Table 1). Declared accuracy of the calculation calls for fur- 
ther consideration of the two-pion final states in the reaction 7* + iV — > iV + tt + 71 (5 
line in Table 1). Indeed, the comparison of the total photoabsorption cross section, de- 
rived from the isobar model, with experimental data shows that in the resonance region at 
2 GeV > W > 1.5 GeV theoretical photoabsorption cross section is less than experimental 
one by the value 100 fib approximately. Two-pion final states in the resonant reactions of the 
form (20) are taken into account in the construction of the cross sections <Jt.l through the 
use of a two-stage model (21). Corresponding contribution to the correction (10) is equal 
to (-6 fieV). But there exists nonresonant contribution of two-pion final states to the total 
cross section <Jtot{Q 2 , W). We included in Table 1 approximate estimate of this contribution 
equal to (-6 fieV). It is based on the assumption that in the whole range of variables Q 2 , 
W: < Q 2 < 1, 1.5 < W < 2 determining the value of the correction (10) the total cross 
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TABLE I: Proton polarizability correction in the Lamb shift of electronic and muonic hydrogen. 



Contribution of the reaction (13) 
to the correction 


(e~p + ) Hz 


(fj, p + ) /ieV 


15 


25 


15 


25 


Contribution of iV7r-states 


- 86.81 


-10.851 


-103.3 


-12.91 


Contribution of iV^-states 


-0.02 


-0.003 


-0.6 


-0.08 


Contribution of .fT-mesons 


-0.03 


-0.004 


-1.1 


-0.14 


Contribution of iW7r-states 


- 0.54 


-0.068 


-12.0 


-1.50 


Nonresonant contribution 


-0.44 


-0.055 


-12.0 


-1.50 


Contribution of subtraction 
term 


0.99 


0.124 


18.4 


2.30 


Summary contribution 


-86.85 


-10.86 


-110.6 


-13.83 



section of virtual photoabsorption is less than experimental data by the value of order 100 
fieV as for the cross section at Q 2 = 0. 

There exists a number of theoretical uncertainties connected with quantities entering 
in the correction (10). In the improved isobar model 0, Hit containing 14 resonances, 
we can omit theoretical error which arises due to the insertion of other high-lying nucleon 
resonances. On our sight the main theoretical error is closely related with the calculation 
of the helicity amplitudes Ai/ 2 (Q 2 ), A 3 / 2 (Q 2 ), Si/ 2 {Q 2 ) in the quark model based on the 
oscillator potential j46j. Only systematical experimental data for the helicity amplitudes 
of the photoproduction on the nucleons Ai/ 2 {0), A 3 / 2 (0) are known with sufficiently high 
accuracy to the present 24]|. In the case of amplitudes for the electroproduction of the 



nucleon resonances experimental data contain only their values at several points Q 2 . So, 
we have no consistent check for the predictions of the oscillator model. Possible theoretical 
uncertainty connected with the calculation of amplitudes Ai/ 2 (Q 2 ), A 3 / 2 (Q 2 ), Si/ 2 (Q 2 ) with 
the account of relativistic corrections can attain the value of order 10 %. Then the theoretical 
error for the correction (10) in the resonance region comprises 20 % from the obtained value, 
that is ±2 fieV for the energy level 2S in muonic hydrogen. There is theoretical uncertainty 
in the contribution due to two-pion nonresonant processes which is presented above. The 
error in this case can constitute no less than 30 % from corresponding contribution that is 
±2 fieV. The essential part of the theoretical error is connected with the subtraction term 
in the dispersion integral (4). Indeed, the increase of the world average value for the proton 
magnetic polarizability from 1.6 x 10~ 4 fm 3 to 1.9 x 10~ 4 fm 3 [24| during last years leads to 
the decrease of the summary contribution to the Lamb shift (2P — 2S) by 0.4 [ieV . The error 
of Pm indicated in Ref.jlij, gives the uncertainty ±0.6 fieV to the theoretical result for the 
shift AEp^(2P — 2S). The obtained value of the proton polarizability contribution to the 
Lamb shift (2P — 2S) in muonic hydrogen is equal to (13.8 ± 2.9) U§V- It is intermediate in 
the value between the result (16 4- 17) fieV calculated in Refs. lg, 17] and 12 ueV, obtained 
in Ref . [l8l| . Our theoretical uncertainty is slightly higher than in Refs. 

0,1 1121 [18j], because 

of the presence of a number of additional theoretical quantities in the isobar model with the 
definite theoretical uncertainties. In the case of electronic hydrogen the value of the proton 
polarizability contribution to the Lamb shift is determined for the most part by the process 
7* + iV — > 7i + iV including the resonance contribution. The obtained shift (—87) Hz of 
the IS" energy level in (e~p + ) is also in the agreement with the previously derived results 
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(-72) Hz in Ref.Q, (-71) Hz in Ref.Q and (-95) Hz in Ref.[u|. 
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